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^ : ABSTRACT 

m 

' Galactic winds and mass outflows are observed both in nearby starburst galax- 

' ies and in high-redshift star-forming galaxies. In this paper we develop a sim- 

, pie analytic model to understand the observed superwind phenomenon. Our 

model is built upon the model of McKee & Ostriker (1977) for the inter- 
stellar medium. It allows one to predict how properties of a superwind, such 



o 

' as wind velocity and mass outflow rate, are related to properties of its star- 

^ \ forming host galaxy, such as size, gas density and star formation rate. The 

model predicts a threshold of star formation rate density for the generation 
' of observable galactic winds. Galaxies with more concentrated star formation 

' produce superwinds with higher velocities. The predicted mass outflow rates 

are comparable to (or slightly larger than) the corresponding star formation 
rates. We apply our model to both local starburst galaxies and high-redshift 
Lyman break galaxies, and find its predictions to be in good agreement with 
current observations. Our model is simple, and so can be easily incorporated 
into numerical simulations and semi-analytical models of galaxy formation. 
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1 INTRODUCTION 

Galactic-scale bulk motions of gas, such as galactic winds and mass outflows, have been ob- 
served both in local starburst galaxies (Forbes et al 2000; Heckman et al 2000 and references 
therein) and in high redshift galaxies (Pettini et al 2000, 2001; Dawson et al 2002; Adelberger 
et al 2002). Evidence for such superwinds and mass outflows comes from X-ray emission pro- 
duced by the hot gas or optical line emission from the warm gas associated with star-forming 
galaxies (Martin, Kobulnicky & Heckman 2002; Strickland et al 2002). Based on medium- 
resolution spectra of Na D (AA5890, 5896) and X-ray observations of some local starburst 
galaxies, Heckman et al (2000) constrained wind velocities in the range ~ 400 — SOOkms"^ 
and mass outflow rates comparable to star formation rates (hereafter SFRs) . Using observa- 
tions of nebular absorption lines and Lj-a emission lines, Pettini et al (2000) and Adelberger 
et al (2002) showed that galactic winds are common in high-redshift Lyman break galaxies. 
The inferred wind velocities for these objects range from several hundred kms~^ to about 
1000 km s~^, with a median value of about 500kms~^. For the bright gravitationally lensed 
Lyman-break galaxy, MS 1512-cB52 {z = 2.72), Pettini et al (2000, 2001) estimated the 
mass outflow rate to be ~ 6OM0yr~^, comparable to the star formation rate of this system 
derived from its UV luminosity. 

These bulk motions of gas (superwinds) in star-forming galaxies are believed to be driven 
by the kinetic energy from supernova (hereafter SN) explosions, and may change the thermal 
and chemical properties of the intergalactic medium (IGM) that forms galaxies (Ferrara, 
Pettini & Schekinov 2000; Furlanetto & Loeb 2002). It is therefore essential to have a proper 
understanding of these phenomena in order to understand galaxy formation itself. 

So far, our understanding of galactic superwinds is quite incomplete. Based on the SN 
remnant evolution model of McKcc & Ostriker (1977), Efstathiou (2000) established a mul- 
tiphase interstellar medium (ISM) model in which galactic winds result from the expansion 
of hot phase gas. Combined with other physical prescriptions such as the infall of the cooling 
gas and star formation, he investigated the formation and evolution of disk galaxies in this 
model, and studied the dependence of mass outflow from a galaxy on the circular velocity 
of its host halo. Silk (2001, 2002) proposed an analytic model for the feedback process and 
outflow in galaxies based on disk gravitational instability and multiphase ISM. Numerically 
there have also been attempts to use simulations to understand how galactic winds can be 
produced (e.g. Tomisaka & Bregman 1993; Suchkov et al 1994; Mac Low & Ferrara 1999). 
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However, current simulations still lack the resolution to deal with processes such as the 
evolution of supernova remnants and the formation of multi-phase gas that may be essen- 
tial for the formation of superwinds. Because of this, simple assumptions are usually made 
about superwinds when studying their impact on the IGM (e.g. Strickland & Stevens 2000; 
Scannapieco, Ferrara & Broadhurst 2000; Scannapieco, Ferrara & Madau 2002; Theuns et 
al 2002; White, Hernquist k Springel 2002). 

In this paper, we attempt to construct a simple analytic model for superwinds. Our model 
is based on the model of McKee & Ostriker (1977) for supernova evolution in the ISM. Our 
goal is to understand how properties of superwinds, such as wind velocity and mass outflow 
rate, are related to the properties of star-forming galaxies, such as size, gas density and SFR. 
The establishment of such relations will allow one to include superwinds in simulations and 
semi-analytical models of galaxy formation. We test our model by comparing its predictions 
with observations of superwinds in local starburst galaxies and in high-redshift Lyman-break 
galaxies. 

Our basic model for the generation of superwinds is similar to that of Efstathiou (2000), 
as both his paper and ours are based on McKee & Ostriker (1977). While Efstathiou focused 
on the details about the assembly and star formation in two example galaxies (one dwarf- 
type and the other Milky- Way type) , we focus on the statistical properties of superwinds in 
star forming galaxies based on simple empirical models about star formation. In addition, 
we compare our theoretical predictions with observations for local starburst galaxies and 
high-redshift Lyman-break galaxies. 

2 MODEL 

In this section, we will use the model of McKee & Ostriker (1977, hereafter M077) for 
supernova evolution in the interstellar medium to construct a model of galactic winds and 
mass outflows for star forming galaxies. Before going into the details, we summarize here 
the basic idea in our modehng (see also Efstathiou 2000, hereafter E2000, for more details). 

In the model of M077, cold star-forming gas in the ISM is assumed to be in the form of 
cold clouds surrounded by the warm ISM with temperature of about ^ 10*^K. Massive stars 
evolve and explode in a few million years as supernovae (SNe), and SN remnants propagate 
to form a low-density hot medium with a temperature of about 10^~^K. The expansion of the 
SN remnants compresses and sweeps out material from the cold clouds in the shock front to 
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produce large-scale bulk motions in the ISM, which can be accelerated by the pressure in the 
SN remnants. We assume such bulk motion to be responsible for the observed superwinds 
and mass outflows. In such a scenario, the properties of the superwind from a galaxy are 
expected to depend not only on the properties of the ISM, but also on other properties of 
the galaxy in consideration. In what follows, we quantify the main processes involved in our 
model. 

2.1 Supernovae evolution, interstellar media and galactic winds 

Star formation generally takes place in giant molecular clouds which consist of many small 
and dense clouds. The differential number density distribution of clouds as a function of 
cloud radius is assumed to follow a power-law (M077, E2000) 

dNc -4 ^ ^ /.N 

——^noa , oi < a < au, (1) 
da 

where a\ and are lower and upper limits for the cloud radii, and Uq is a normalization 
constant. We take, following M077, 

ai = 0.5pc, Ou/oi = 20. (2) 

We caution that the exact value of the lower limit radius (ai) is unclear. It is generally 
believed to be in the range of 0.5 to Ipc (M077); recent observations by Olmi & Testi 
(2002) suggested that a\ ~ 0.55 pc. 

It is easy to show by integration that the local cloud number density, rid, is related to 
no by 

^ci = '«'o/3af, for 3> oi. (3) 

Assuming that Pc is the mean mass density within individual clouds, the mean cold gas 

density pcoid can be written as 

dNc 47r o _ , 47r _ , , ^ , 

Pcoid = / —J —a Pcda = —pcTioln—. (4) 

Jai da 6 6 ai 

Thus, the volume flUing factor of the cold gas is 



J. Pcold / f. \ 

Jc = 5) 

Pc 

The porosity P is related to /c by fc = e~^. 

According to M077, an expanding SN remnant will sweep up the ambient cold gas with 
a total mass of 

Mev - 540 Esi'/^nh'/'E-^/^Me, (6) 
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where E51 is the energy output by one SN in unit of lO^^erg, rih is the hot gas density 
interior to the SN remnant and E is the evaporation parameter 

^ ^ 7 _ 'jaiHaJai) 

47rainci0fc (pk Pcoid ' 

Here 7 is the ratio of the blast wave velocity to the isothermal sound speed of the hot phase, 

which is usually taken to be 2.5 (E2000), and the parameter 0^ denotes the efficiency of 

conduction relative to the classic thermal conductivity of the clouds, (pk is roughly in the 

range of 0.1 to 0.01 due to the suppression of heat conduction by tangled magnetic fields, 

turbulence, etc (E2000). 

As argued by several authors (Silk 1997, 2001, 2002; E2000; Clarke & Oey 2002), star 
formation in galaxies may proceed in such a way that the porosity P is always maintained 
as a constant value close to unity in star formation regions. This constancy in the porosity P 
is supported by observations in the Milky Way, and can be physically understood as follows. 
If the porosity were too large, the cold gas fraction would be small and the star formation 
activity would be reduced. The hot gas will then cool quickly, and hence leading to a lower 
porosity. Conversely, if the porosity were too small, the cold gas fraction would be large. As 
the dynamical timescale is likely to be short in the star formation regions, star formation 
will ensue and hence lead to a higher SFR and higher porosity. In the present study, we 
assume the porosity to be unity everywhere in a star formation region. 

Adopting Ou/oi = 20, we can rewrite the evaporation parameter (E) as 

Notice that the value of E is not very sensitive to the ratio Ou/oi because it enters eq. (8) 
logarithmically. For convenience, we define a new parameter 

= Eo = 95pc2 (9) 

which is the evaporation parameter normalized to the value close to the solar neighborhood 
(E2000). Numerically, we find 

Following M077 and E2000, we can also obtain the temperature of the hot gas 
Th = 6.6 X 105(5_i3 E5i/s/7)°'' K, (11) 
and the mass evaporation rate per unit volume 

Pev = 2.7 X 10-''S'4h'-'' E5i°-^VE°-''Me pc" (12) 
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respectively, where S'_i3 is defined as the the SN explosion rate in unit of 10~^^pc~^yr~^. 
From equations (11), (12) and the definition of /s [eq. (9)], we see that the temperature Th 
and the mass evaporation rate pev are in fact independent of the parameter 7. Note that 
Pev increases with supernova rate density (which is proportional to the SFR density) as a 
power of 0.71. Thus, for a given total star formation, the total outflow rate increases with 
the volume of the star formation region. 

The hot phase will first expand to form super bubbles which compress the ambient ISM. 
The compressed ambient ISM will be driven to form a "wind" with a velocity comparable 
to the isothermal sound speed of the hot phase, which is given by 

d = (A;Th//xmp)^/2 = 37T^^^ kms'^ (13) 

with T5 the hot gas temperature in unit of lO^K and fj, the mean molecular weight per 
particle; we take /i — 0.61, a value appropriate for a fully ionized primordial gas. Because of 
the conservation of the specific enthalpy, the wind will accelerate and reach a bulk terminal 
speed to form a superwind (E2000). The terminal wind velocity is related to the isothermal 

sound speed by 

^'wind = ^wCi- (14) 

We take ^ V^23, instead of a value -\/5; our value is appropriate when part of the thermal 
energy is lost radiatively [see Appendix B in E2000 for details]. Observationally, only when 
the initial sound speed of the hot gas is larger than ~ lOOkms"^, can galactic winds be 
readily observed (e.g. Heckman et al 2000). From eq. (14), this sound speed corresponds to 
roughly 160kms~^ in the terminal velocity. 

Note that dark matter halos, which dominate the potential wells of galaxies, have not 
entered our discussions explicitly. Thus, the local properties of outfiows, such as wind velocity 
and evaporation rate per unit volume, depend only on the properties of the ISM in the star- 
forming region. This is consistent with the observations of X-rays from hot gas in galaxies 
and the detailed analysis of wind velocities by Heckman et al (2000), who found that terminal 
velocities of the superwinds are almost independent of the galaxy potential wells. However, 
the dark halo of a galaxy will play an important role in determining whether the outfiow 
can eventually escape from the galaxy or fall back into the galaxy. 

It is clear from equations (11) and (12) that the model predictions for the wind properties 
(velocity and mass outfiow rate) depend on the choices of the values of 0fc and ai. As discussed 
above, the values of these two parameters are probably in the ranges 0.01 ;^ 0jk ^ 0.1 and 
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0.5 pc ^ a\ ^ 1 pc. Using equations (11) and (12), one can estimate that these ranges 
correspond to a factor of ^ 2 in the predicted wind velocity, and a factor of ;^ 3 in the 
predicted mass outflow rate. These factors should serve as an indication of the uncertainties 
in our model predictions. 

One assumption we made is that the porosity P retains close to unity. While this seems 
to be supported by observation but hard to prove rigorously (Silk 2001), it is important 
to examine the dependence of our resuhs on the value of P. In general, galactic winds are 
more efficient during the early adiabatic stage of SN remnant evolution, when the porosity 
P ^ 1. In this case, SN remnants do not overlap with each other. The physical prescriptions 
we developed for P = 1 are valid for P < 1, except that the wind velocity I'wind ^-nd 
mass outflow rate Pev must be modified by a factor of P~'^/'^ and (1 — e~^)/[P^/^(l — e~^)], 
respectively (M077). For a range of P from 0.2 to 2, the wind velocity changes by less than 
25% while the mass outflow rate changes by less than 20% from our canonical results (for 
P = 1). Clearly both Wwind and pev are insensitive to P. As SN remnants evolve, they may 
overlap with each other and P may becomes larger than 1. In this case, the model will break 
down. As discussed above, in real galaxies, P may be regulated so that its value cannot 
be much larger than 1. We therefore beheve that the assumption P = 1 does not lead to 
significant error in our results. 

2.2 Star Formation in Spherical Regions 

So far we have only considered superwinds and mass outflows launched locally from a particu- 
lar location in a star formation region. The derived superwinds and mass outflows therefore 
only depend on the local star formation activities. The observed wind velocity and mass 
outflow rate from a galaxy should be a proper average and sum of the wind velocity and 
outflow rate over the whole star formation region. It is these global quantities that can be 
directly compared with observation. In this subsection, we consider a simple model where 
star formation takes place uniformly in a spherical region. Here the wind and outflow are 
simply evaluated by the SFR density and the integration over the volume within star form- 
ing region. We treat individual star forming cells as source of winds but ignore the detailed 
transfer of energy and mass from the center to the outer region. Our treatment is obviously 
simplistic; a more careful consideration must consider wind interactions which result in the 
eventual outflow from the surface of the spherical region. 
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For a spherical region of uniform star formation, S'_i3 can be written as 

.2xio»f_* uj^ym", (15) 



where is the SFR within radius R and Mpg is the stellar mass of star formation per 
supernova explosion. For a Salpeter mass function {n{M)dM oc M~^-^^o?M), Mpg ^ 125M0 
with a stellar mass range from 0.1 to SOMq. 

Using equation (14), we can estimate the galactic wind terminal velocity 



Vwind = 623 km s ^ 



lOOMeyr-V VI25M0/ V^kpc^ 



IC 



■0.29 



(16) 



and the corresponding total mass evaporate rate 

= 133M,y..-' I r fiV°" f^y\-o-. (17) 

\^lOOM0yr-i; yUhMQ) ^Ikpc/ ' ^ ' 

For convenience we have defined /C as a function of /s and 7 as 

K=f^V^V', (18) 



^21.5 J V2.5, 

which gives the dependence of I'wind a-nd Mey on the properties of the ISM. 

For fixed /C, the wind velocity depends only on the SFR density, M^/R^, because the 
temperature of the hot gas, hence the wind velocity of the driven ISM, depends only on the 
SN explosion rate density [cf. eq. (12)]. The wind velocity increases with the SFR density. 
Thus, galaxies with a more compact cold gas distribution will produce winds with higher 
velocity, no matter what the total amount of cold gas and total SFR are. This implies 
that superwinds are more likely to be observed in local starbursts and high-redshift star 
forming galaxies where the gas density is high and star formation activity is concentrated. 
For a given SFR density, the wind velocity increases with the lower limit of the cloud radii, 
because larger clouds leads to a larger value of /e [cf. eq.(9)]. The wind velocity also increases 
with decreasing 0^, because a higher 0^ implies a higher energy loss due to thermal transfer. 
As discussed in the last subsection, we adopt the lower limit of the wind velocity to be 
~ 160kms~^ which, according to Heckman (2000), is the lower limit that can lead to an 
observable superwind. This limit directly translates into a lower limit on the SFR density, 
0.01 — 0.1 M0yr~^ kpc~^, below which no observable wind is expected. 

On the other hand, the mass outflow rate has an additional dependence on the size of 
star formation region. In Fig. 1, we show the predicted mass outflow rate as a function 
of SFR for three sets of parameters. The solid, dashed and dotted lines denote the results 
for (0fc,ai) equal to (0.1, 0.5 pc), (0.01, 0.5 pc), and (0.01, Ipc), respectively. The thin and 
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0.01 0.1 1 10 100 1000 

M. (M^ yr-0 

Figure 1. The predicted mass outflow rate as a function of the SFR for three sets of parameters, shown as the soUd, dashed 
and dotted hues corresponding to (t/ife, ai) equal to (0.1,0.5 pc), (0.01, 0.5 pc), and (0.01, 1 pc), respectively. The thin and thick 
lines are for two different sizes of star formation regions, 0.1 kpc and Ikpc, respectively. 

thick lines show resuhs for the two adopted size of star formation region, 0.1 kpc and 1 kpc, 
respectively. For a given total SFR, a larger star formation region gives a larger mass outflow 
rate, because the increase of the outflow rate density with SFR density is slower than linear 
(see eq. 12). Note, however, this increase of M^^ with size cannot go indefinitely, because 
our assumption that P = 1 will fail when the gas density becomes so low that the volume 
is too large for the SN remnants to fill. 

Since the total power in the wind, which is proportional to Me^v^^^^, should be propor- 
tional to the total SFR, the dependence of Mev on model parameters can be understood 
in the same way as the dependence of wind velocity on model parameters. From Fig. 1 we 
see that the predicted mass outflow rates are comparable to (or higher than) the corre- 
sponding star formation rates for a wide range of model parameters. This is consistent with 
observation, e.g. Heckman et al. (2000) and Pettini et al. (2000; 2001; 2002). 

2.3 Star Formation in Galactic Disks 

In this subsection, we consider another case where star formation is assumed to be in thin 
galactic disks. This may be more realistic than the spherical model described above, because 
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the cold gas supporting star formation generally has angular momentum and is found in a 
rotation-supported thin disk. 

For simplicity, we assume the velocity dispersion (cTg) of the cold clouds in a disk to be 
constant everywhere. Throughout this paper, we choose Ug to be 10 km s"^, a value consistent 
with observation (Stark & Brand 1989). The equation of hydrostatic equilibrium gives the 
following solution of cold gas density distribution in the vertical direction: 

= ^r"" (A) • 

where /ig is the cold gas surface density and Hg is the scale height given by 

= ^ = 74 f ,^r^V' fr^V PC, (20) 

^ vrG/Xg V^epc-^y Vkms-V ^' ^ ' 

and G is the gravitational constant (Spitzer 1942). 

Kennicutt (1998) studied star formation in a wide range of physical conditions, ranging 
from quiescent galactic disks to starburst regions. He derived an empirical law for the SFR 
per unit area as a function of the cold gas surface density, 

A. = 2.5 X 10-^° (77^1 ' ' Meyr-^ pc'^. (21) 
It then follows that the SFR density is 



2.4 



= 1-7 X 10-" ^-^) Meyr-'pc- (22) 



-2 



mg V^opc-^ Vkms-1 

and the SN explosion rate S-vi is 

^-13 = 10^^ X = 1.35 X 10- [ -^"l " {-^\ " [ -^"l " , (23) 
Mps VM0pc-2y Vkms-V \i25M0y ' ^ ' 

where Mpg is again the mass of stars formation corresponding to one supernova explosion. 

Substituting the above equation into eqs. (11) and (14), we can infer the wind velocity 
f wind to be 

The corresponding mass outflow rate per unit area can be obtained by multiplying the mass 
outflow rate per unit volume and the scale height 

where /C is deflned in eq. (18). 

As one can see, for given Ug and /C, the wind velocity increases with \x^. Thus, superwinds 
are more likely to occur in compact systems of cold gas. The wind velocity decreases with 
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Figure 2. The predicted mass outflow rate per unit area as a function of the SFR per unit area. The solid and dashed hues 
are the results where the lower limit of the cloud radius is chosen to be 0.5 pc and Ipc, respectively. The thin and thick lines 
denote results with the parameter (jik chosen to be 0.01 and 0.1, respectively. 



ag, because a higher value of ag imphes a larger disk thickness, and hence a lower density of 
cold gas. The mass outflow rate per unit area also increases with gas surface density, because 
of the dependence of the SFR per unit area on the surface density of cold gas. 

As for the spherical star formation region in last subsection, we adopt the same threshold 
for the wind velocity, below which winds may not yield observable signatures (Heckman et 
al. 2000). This threshold corresponds to a lower limit on the SFR per unit area, which can 
be obtained from eq. (24). For the reasonable ranges of (pk and ai, i.e., (pk = 0.01 — 0.1 and 
ai = 0.5 — 1 pc, the limit is about 



which, from eq. (21), corresponds to a cold gas surface density about 10 — 40MqPc~^. This 
limit is consistent with the observational estimation quoted by Heckman (2001). 

The predicted mass outflow rate per unit function of the star formation rate 

per unit area is shown in Fig. 2. The solid and dashed lines are results when the lower limit 
of the cloud radius is chosen to be 0.5 pc and Ipc, respectively. The thin and thick lines 
denote results with the parameter 0^ chosen to be 0.01 and 0.1, respectively. For any given 
(pk and ai, the mass outflow rate per unit area increases with the SFR per unit area, as 
pointed out before. Note that the predicted mass outflow rates per unit area are comparable 
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to (or higher than) the SFRs per unit area, which is consistent with observation (Heckman 
et al 2000; Pettini et al 2000; 2001; 2002). 

So far we have considered only local properties of winds in galactic disks. To facilitate 
direct comparison between model predictions and observations, we must consider globally- 
averaged quantities of the wind. As a simple example, we assume that the cold gas is dis- 
tributed in an exponential disk 

lig = //oexp(-r/rd), /Iq = 77-^ (27) 

where /xq is the central surface density of cold gas, Mg is the total cold gas mass and is 
the scale length for the disk. Note that our results are not sensitive to the assumed cold gas 
surface density distribution because the cold gas surface density enters the wind velocity 
and mass outflow rate only with moderate power-law indices (0.35 and 0.7) in eqs. (24) and 
(25), respectively. 

As discussed above, there is a threshold in the SFR per unit area (eq. 26), i.e., the lower 
limit of the cold gas surface density, above which observable galactic winds can be produced. 

If the central surface density /xq for a disk galaxy is lower than this, no observable galactic 
wind will occur. When /io is above the threshold, the threshold implies a critical radius rcr 
for any given fC, beyond which no outflow will occur. The critical radius satisfies 

A*cr = A*oexp(-rci./rd). (28) 
It is easy to calculate the total mass outflow rate within this critical radius 

= 6.2X10-M,yr.(ii^)^(^)°' 

Vkms-V Vl25M0y ^ '^"^ 

where the function J^{x, /Iq) is defined as 

T{x,iio) = 1 - 



1 — xln ( 



!F{x, fjLo) describes the fraction of area that contributes to the observable galactic wind. For 
any given x, T will first increase with /xq for small values of //q, and then decrease with 
increasing /xq for large values of /Xq. 
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Figure 3. The predicted mass outflow rate weighted wind velocity as a function of /ig where the soHd and dashed lines have 
the same notations as that in Fig. 2 with the dotted horizon line marking the lower limit of the velocity, ~ 160kms~^. Here 
the velocity dispersion of the cold gas clouds is assumed to be lOkms"^. 

The observed superwind velocity is an appropriate average of the wind velocities over 
the whole disk. Since the mass outflow rate per unit area is proportional to /Xg^, we define 
the global wind velocity to be the average of the wind velocities weighted by /ig^. This wind 
velocity is 

- ^ {^r i^r m°'^^ (^) '^^^ 

In reality, the observed wind velocity depends on where the wind originates, and our weight- 
ing scheme can only serve as an approximation. We have also made calculations using a 
weight proportional to fig and to /Xg^ (i.e. to the surface density of SFR). The results for 
these three weighting schemes differ only by 20%. Note that for given JC the wind velocity 
for an exponential disk depends only on the central surface density of cold gas, while the 
total mass outflow rate for the disk depends in addition on the scale length (or the total 
cold gas mass). 

Fig. 3 shows the the predicted global wind velocity for an exponential disk with fiQ ^ ficr 
as a function of cold gas central surface density /io for four different combinations of (pk 
and ai. The notation of the figure is the same as that in Fig. 2 with the dotted horizon 
line marking the lower limit of the velocity, ~ 160kms~^, as discussed above. For given 0^ 
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Figure 4. The predicted mass outflow rate (solid lines) and the corresponding star formation rate (dashed lines) as a function 
of ^0 for given the total cold gas mass, where lines from top to bottom denote the total cold gas mass to be 2 X 10^^ Mq, 
2 X IO^^Mq and 2 X IO^Mq respectively. Here the velocity dispersion of the cold gas clouds is assumed to be lOkms"^, and 
(a) and (b) denote the results of (cjifc, ai) equal to (0.01, 1 pc) and (0.1, 0.5 pc), respectively. 

and a\, the predicted global wind velocities increase with the cold gas surface density /xq- 
Galaxies with low cold gas surface density (e.g., low surface brightness galaxies or high- 
surface brightness galaxies with low cold gas surface density such as our Milky Way) will 
not produce observable galactic winds because their star formation rate is too low. 

As an example, Fig. 4 plots the predicted total mass outflow rate (solid lines) and the 
corresponding SFR (dashed lines) as a function of /io with (a) and (b) denoting the results of 
(0fc, CLi) equal to (0.01, Ipc) and (0.1, 0.5 pc), respectively. The lines from top to bottom are 
results for a total cold gas mass to be 2 x IO^-'^Mq, 2 x 10^°Mq and 2 x IO^Mq, respectively. 
From the figure we see that the cold gas central density /io must be larger than a critical 
value /icr (about 10 and 4OM0 pc~^ for the values of 0^ and ai adopted here respectively) 
for an observable wind to occur. If /iq < /icr, no observable galactic wind will be produced 
no matter how large the total cold gas mass is. 

For /io ^ /icr, the total mass outfiow rates increase with the total mass of cold gas for a 
given /io because of the increase of the SFR and the star formation region which can produce 
observable outfiows. For a given total cold gas mass, the predicted mass outfiow rate first 
increases rapidly with /io, because the gas mass that can produce observable winds [i.e. rcr 
in eq. (28)] increases with /io. The mass outfiow rate as a function of /io reaches a plateau 
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and then decreases when /iq is larger than ~ IO^MqPc"^. This happens because when /io 
is high enough, the total gas mass that can produce observable galactic winds is saturated 
while the SFR density increases with increasing /xq- This increase in SFR density reduces 
the mass outflow rate, again because the increase of Pev with SFR density is slower than 
linear (see eq. 12). Note that the predicted mass outflow rates are comparable to the SFR 
for a very wide range of iiq, which is consistent with current observations. 



3 COMPARISONS WITH OBSERVATIONS 

As discussed in the introduction, many observational studies have investigated galactic winds 
both in local starburst galaxies and in high-redshift star forming galaxies. In this section, 
we examine whether the predictions of our model can match the current observational data. 

For definiteness, we will adopt 0^ = 0.01 and ai — 1 pc throughout this section. As we 
discussed in Sec. 2.1, other plausible choices of these parameters may result in a lower wind 
velocity by a factor of ^ 2 and a higher mass outflow rate by a factor of ^ 3. The choice is 
quite arbitrary, since the exact values of 4>k and ai appropriate for star forming galaxies are 
not known a priori. Our chosen values do give reasonable agreement with the observational 
results to be discussed below. 

3.1 Local starburst galaxies 

Based on the Na D absorption hues, Heckman et al (2000) estimated galactic wind velocities 

for local starburst galaxies. They found that the wind velocities are in the range of 400 to 

800kms~^. They also estimated mass outflow rates, which they claimed to be comparable 

to the SFRs. The observational results for the mass outflow rates are only qualitative and 

we have shown that they are consistent with our model predictions. In what follows, we will 

primarily focus on the wind velocities which are more accurately determined. 

At first we assume that the star formation takes place within spherical regions. We can 

then estimate the global galactic winds for these galaxies using equation (16), provided that 

their SFRs and sizes are known. Wc estimate the SFRs for these starburst galaxies using 

their FIR luminosities (Kennicutt 1998) 

SFR ^ -^FiR /„2^ 
M0yr-i 5.8x109Lq' ^ ' 

The sizes of the star forming regions for the galaxies in consideration are not available. As 
a rough estimate, we take the sizes to be the same as the observational slit widths, which, 
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Figure 5. The predicted wind velocities vs observed values for local starburst galaxies. The crosses denote the observations of 
Na D absorption lines and open squares denote the estimated results based on the the hot gas temperatures from ROSAT and 
ASCA. The lines denote the agreement between observation and model predictions, (a) for spherical star formation region; (b) 
for exponential star formation disk (see the text for details). 

according to Heckman et al (2000), is a reasonable match to the typical sizes of powerful 
starbursts. The observational data are from Table 4 in Heckman et al (2000). The predicted 
wind velocities vs the observed ones are plotted as crosses in Fig. 5(a). The lines in the 
figure denote the agreement between observation and model predictions. As can be seen, 
our model predictions roughly match the observations within the model and observational 
uncertainties. 

Another way to empirically estimate the super wind velocities for galaxies is based on 
X-ray observations. Galactic winds will spread out the hot gas within galactic halos and 
contribute to the soft component of extended X-ray emission. As the hot gas temperature 
can be inferred from X-ray observations (e.g., from ROSAT or ASCA), we can estimate the 
wind velocities from the temperature of the soft X-ray components using equation (11). The 
galaxies we selected are similar to those in Heckman et al (2000). They are NGC4449 (Delia 
Ceca, Griffiths & Heckman 1997), NGC2146 (Delia ceca et al. 1999), NGC253, NGC3079, 
M82, NGC4631 (Dahlem, Weaver & Heckman 1998), NGC1569 (Delia ceca et al. 1996), 
Arp299 (Heckman et al 1999), NGC6240 (Iwasawa & Comastri 1998), NGC1808 (Awaki et 
al. 1996). Note that we assume the relation between wind velocity and soft X-ray temperature 
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to be Kvind = \/'2.-^kTx rather than y/5kTx, the relation adopted by Heckman et al (2000). 
The reason for our choice is discussed below eq. (14). 

We can estimate the SFRs of these galaxies based on their FIR luminosities according to 
equation (32). The sizes of their star formation regions are adopted from the corresponding 
Ha observations, which can be found in the references listed above. The predicted wind 
velocities can then be obtained using eq. (16). The open squares in Fig. 5(a) indicate the 
predicted wind velocities and the estimated values from the X-ray observations. We find 
that the predictions match the estimated values reasonably well. 

As pointed out in Sec. 2.3, the predicted results are not sensitive to the assumed gas 
surface distribution if we assume starbursts take place in exponential disks. Assuming that 
the cold gas in these galaxies is distributed exponentially and that observed star formation 
region contain just one half of the total star formation activity in the galaxies, we can also 
predict the superwind velocities for the above starburst galaxies based on equation (31). 
The model predictions versus observational results are shown in Fig. 5(b). Here again, the 
model predictions match observation. 

3.2 Lyman brezLk galcLxies 

The UV dropout method has been very successful in identifying active star forming galaxies 
(Lyman break galaxies, hereafter LBGs) at a redshift of z ~ 3 (Steidel, Pettini & Hamilton 
1995). Much observational and theoretical work has investigated their physical nature (see 
Mo & Fukugita 1996; Mo, Mao & White 1999; Kauffmann et al. 1999; Katz, Hernquist & 
Weinberg 1999; Shu 2000; Shu, Mao & Mo 2001). Based on observations of Lyman alpha and 
nebular emission lines, LBGs are also inferred to display galactic winds (Pettini et al 2001; 
2002; Adclbcrgcr et al 2002). LBGs display a wide range of wind velocities, 250 to nearly 
1000 kms~^ with a median value about 500-600 km s^^. As for local starburst galaxies, the 
mass outflow rates for LBGs are difficult to establish observationally. Pettini et al (2000; 
2002) measured the mass outflow rate for a specific LBG, the lensed and magnified MS 1512- 
cB58; the value is ~ 6OM0yr~^, comparable to its SFR obtained from its UV luminosity. 

At present it is unclear whether the star formation in LBGs occurs in quasi-spherical 
regions or more in a disk-like geometry. For this reason, we will consider both scenarios 
below in turn. 

If the star formation activities in LBGs are quasi-spherical and homogeneous, then we 
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Figure 6. The predicted distribution of the galactic winds for LBGs. The soUd and dashed histograms denote results where 
star formation is assumed to take place in spherical regions and thin exponential disks, respectively. 

can use equations equations (16) and (17) to estimate the mass outflow rates and galactic 
wind velocities for LBGs if we know the star formation rates and the sizes of their star 
formation regions. Following the procedure outlined in Shu, Mao & Mao (2001), we adopt 
the SFR distribution of LBGs from their dust-corrected UV luminosity and a log-normal 
size distribution. We use Monte Carlo simulations to select galaxies following these SFR and 
size distributions. We then substitute these two quantities into eqs. (16-17) to obtain the 
predicted distributions of the galactic wind velocities and the corresponding mass outflow 
rates. The results are shown as the solid histograms in Figures 6 and 7, respectively. The 
predicted winds have velocities in the range from 200 to 900kms~^ and a median value of 
450 km s^^, in agreement with the quite limited observational results. The predicted median 
value of the mass outflow rate is about 8OM0yr~^, similar to the inferred value for MS 
1512-CB58. 

The second scenario we consider for LBGs is that their star formation occurs in a pre- 
assembled disk. For simplicity, we assume that the cold gas surface density has an exponential 
distribution. In this case, the model discussed in Sec. 2.3 can be applied to the LBG pop- 
ulation. We use a procedure identical to that used in the spherical model to simulate the 
LBG population. 

Based on eqs. (29-31), the estimated outflow rates and wind velocities can be obtained. 
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Figure 7. The predicted distribution of the mass outflow rates for LBGs. The soUd and dashed histograms denote results 
where star formation is assumed to occur in spherical regions and thin exponential disks, respectively. 

The dashed histograms in Figs. 6 and 7 show the model predictions for the distributions of the 
galactic wind velocities and the corresponding mass outflow rates for LBGs, respectively. It 
can be seen that the median value of the predicted wind velocities for LBGs is ~ 600 km s^^. 
The median value of the predicted mass outflow rates is about 60 M0yr^^. Both are consistent 
with observations. 

4 ESCAPING GALACTIC WINDS AS A FUNCTION OF Vc AND 
REDSHIFT 

As we discussed in Sec. 2, galaxies with more compact cold gas distributions (and hence 
more compact star formation activity) will produce stronger galactic winds. Galactic winds 
will thus occur preferentially in local starburst galaxies and high redshift galaxies. If the 
speed of a galactic wind is smaller than the escape velocity of the host galaxy, the mass 
outflow will fall back into the galaxy and form a galactic fountain. Otherwise, the outflow 
will escape from the host. This we call an "escaping" outflow. The mechanical energy of the 
outflow can heat the the intergalactic medium (IGM) while heavy elements contained in the 
wind can chemically enrich this gas. It is therefore important and interesting to estimate the 
fraction of galaxies (by number) that will have escaping outflows at different redshifts. We 
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Figure 8. The predicted fraction of galaxies wliere tlie superwinds escape from tlieir dark matter lialocs is shown as a function 
of circular velocity at six different redshifts. In each panel, the four solid lines correspond to a gas faction of mg =0.1, 0.05, 
0.025 and 0.01 (from top to bottom), respectively. The predicted mass fraction F(Vc) in halos for a given Vc is shown as a 
dashed line at each redshift in the current ACDM cosmogony. 

address this question in the currently preferred ACDM cosmogony with a matter density 
fim = 0.3, a cosmological constant Vlx = 0.7, a Hubble constant h = 70kms~^ Mpc~^; the 
power-spectrum is described by a shape parameter F = 0.2 and the usual normalisation 
constant as = 0.9. 

We assume that the distribution of cold gas within a galaxy is exponential; our results 
will not change significantly if we adopt other profiles (see Sec. 2.3). We use Mo, Mao & 
White (1998) to model the formation of disks in dark matter haloes. In this model, the disk 
properties are determined by the halo circular velocity, Vc, a dimensionless spin parameter. A, 
and the fraction of total mass {mg) that settles into the exponential disk. The spin parameter 
follows a log-normal distribution with a median A = 0.05 and a dispersion ainA = 0.5 (e.g., 
Warren et al 1992; Lemson & Kauffmann 1999). Once these three parameters (K, A, and 
mg) are specified, the disk properties and hence the star formation and outflow properties 
can be evaluated using the formalism developed here. In Fig. 8, we show the number fraction 
of galaxies with escaping winds occuring in the disks as a function of circular velocity at 
six different redshifts. For each redshift, four solid lines are plotted, corresponding to four 
different mg values. Note that, for these results, we have averaged over the spin parameter 
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distributions. In the same figure, we also show the mass fraction F{Vc) in halos with a given 
circular velocity Vc in logrithmic bin at different redshifts, i.e., 

F(K)dlogVe = ^-^(>^c,^) ^iogV„ (33) 
dlogVc 

where F{> Vc,z) is the mass fraction in haloes with circular velocity larger than Vc as 
predicted by the updated Press- Schechter formahsm (c.f. Mo & White 2002). From the 
dashed line, one can easily estimate the mass fraction of total halos and the typical halo 
circular velocity which will produce escaping mass outflows. 

As expected, for any given rUg, a larger fraction of small halos produces escaping outflows 
at any rcdshift because of their shallower potential wells. Hence, during their evolution, a 
significant fraction of their baryons will be lost due to the galactic winds. It also implies 
that most small halos in the local universe are dark matter dominated which is consistent 
with observation (e.g. Persic, Salucci & Stel 1996). Given V^, galaxies with larger nig have a 
larger probability to produce an escaping outflow at any redshift because of their more active 
star formation. As the redshift increases, more and more galaxies produce escaping outflows 
because of the increase in their SFR densities. For example at z = 3, more than 80% galaxies 
with Vc ^ 200kms~^ will produce mass outflows that escape for rUg ^ 0.05. Hence escaping 
outflows will occur commonly in LBGs because the median value of their circular velocity is 
about 150kms~^ and their nig may be ^ 0.06 (Shu, Mao & Mo 2001). Hence, a signiflcant 
number of galaxies will contribute to the heating of the IGM or ICM. Subsequent galaxy 
formation therefore will occur in preheated intergalactic media, as argued by Mo & Mao 
(2001). Galactic winds therefore should be taken into account for high-rcdshift galaxies. 

As a comparison, it is interesting to study how the star formation rate is partitioned 
in haloes with different circular velocities. Specifically, we study the differential probability 
distribution of star formation rate in logrithmic scales of V^, i.e., 

^'^')'''°«^'- JSFR(v!)dlogv/ '°^^" ^''^ 
where SFR is the star formation rate in galaxies for a given circular velocity Vc averaged 

over all spin parameter distributions. Note that this distribution function is normalised to 

unity for Vc between 30 km/s and 400 km/s. 

We define the accumulative probability distribution function for galaxies with escaping 

outflows, /ro(< Vc), function of V^, 

f . .,._ /30SFResc(Vc)dVc 
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Figure 9. The predicted fraction of galaxies weighted by the SFR as a function of circular velocity at six different redshifts. 
Gas cooling and feedback effect have been taken into account and the currently favoured ACDM cosmogony is adopted (see 
text for details). 



where SFResc is the star formation rate in galaxies with escaping outflows for a given circular 
velocity Vc- 

To obtain these two probability distributions, we have assumed that all galaxies are disk 
galaxies. Furthermore, the gas fraction rrig in the disks is taken to be the fraction of gas that 
eventually cool in their host haloes; we take the cooling function from Sutherland & Dopita 
(1993) assuming the metallicity of 0.01 solar value. The influence of the feedback on rUg is 
also taken into account by adopting 



as suggested by Dekel & Silk (1986) and White & Frenk (1991). Here mgo is the maximum 
baryon fraction that can cool in the haloes to form disks. 

We show the two probability distributions in Figures 9 and 10 for six different red- 
shifts, respectively. We can see from Fig. 9 that the SFR fraction of galaxies peaks at 
Vc ~ 100 - 200kms~^ at redshift z > 1. However, at lower redshift, most of the SFR is 
contributed by more massive galaxies because the feedback at low redshift is not as efficient 
and most baryons within massive halos can cool to form stars. From Fig. 10, we see that 
the contribution of SFR for galaxies with escaping outflows are dominated by galaxies with 



mgo 



(36) 



m„ = 
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Figure 10. The predicted accumulative probability distribution function of galaxies with escaping outflows as a function of 
halo circular velocity at six different redshifts. For each redshift, the solid and dashed lines denote results for rrtgo =0.1 and 
0.05, respectively. Gas cooling and feedback effect have been taken into account and the currently favoured ACDM cosmogony 
is adopted (see text for details). 

Vc ~ 200kms~^ at redshift z ^ 3, while at higher z significant contribution comes from 
systems with ~ 150kms~^. This is because, ahhough at lower redshift the SFR is shifted 
toward systems with higher V^, escaping superwinds are more difficult to produce in these 
systems. At 2; > 3, more than 40% of all stars are formed in systems where escaping winds 
are expected. Note that a change in the value of mgo does not affect the result significantly. 
At high redshift, the number density of high Vc systems goes down rapidly and so not many 
stars can form in these systems. 

5 SUMMARY 

In this paper, we develop a simple analytic model to understand superwinds and mass 
outfiows that are observed in local starburst galaxies and high-redshift star-forming galaxies. 
Our scenario is based on the model of Mckee & Ostriker (1977) for the evolution of SN 
remnants in the ISM, which is similar to Efstathiou (2000). We show that the properties of 
superwinds, such as wind velocity and mass outfiow rate, depend not only on the properties 
of the ISM, but also on the global properties of star forming galaxies. Our main conclusions 
are: 
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• Observable winds are produced only in systems where the density of SFR is higher than 
some threshold. This is consistent with the observational results given by Heckman (2001) 
and Heckman et al (2000), and implies that low surface-brightness galaxies and high surface- 
brightness galaxies with low cold gas surface density cannot produce large-scale superwinds. 

• The velocity of superwinds driven by SN explosions depends only on the SFR density. 
Galaxies in which current star formation is confined in a more compact region can produce 
winds with higher velocity. This implies that superwinds are expected in local starbursts and 
in high-redshift star-forming galaxies. The mass outflow rate depends, in addition, on the 
size of star formation region; a larger star formation region allows more mass to be loaded 
in the wind. 

• The predicted mass outflow rates are comparable to or higher than the corresponding 
SFRs, consistent with current observations both for local starburst galaxies (Heckman et al 
2000) and for high redshift star-forming galaxies (Pettini et al 2000; 2001; 2002). 

• The predicted wind velocity and outflow rate have no explicit dependence on the prop- 
erties of dark halos which dominate the potential wells. Thus, galactic winds and mass 
outflows can occur in a variety of halos, provided that the cold gas density is high enough. 
This is in agreement with the observation that there is no strong correlation between wind 
properties and the total mass of a galaxy (Heckman et al 2000). The potential well of a 
galaxy will, however, determine whether the outflows escape from the galaxy or return. 

• The fraction of galaxies with superwinds that will eventually escape from their dark 
matter haloes is a function of circular velocity and redshift. Our model predicts that the 
fraction is high in low circular velocity systems such as dwarf galaxies although their con- 
tribution to the total SFR is small. More interestingly, we find that the fraction is high for 
galaxies at high redshifts, such as the LBGs at redshift z — 3. The winds will undoubt- 
edly heat and chemically contaminate the IGM, and hence have important implications for 
subsequent galaxy formations in the preheated IGM. 

We apply our model to make predictions for the properties of the winds expected from 
local starburst galaxies and high-redshift Lyman-break galaxies. These predictions can match 
many of the observed properties. We therefore believe that our model catches the main points 
required to model the superwind phenomenon. Our model is also simple, and so can be easily 
incorporated into numerical simulations and semi- analytical models of galaxy formation. 
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